Role of Chiral Symmetry in the 
Nuclear Many-Body Problem^ 



N. Kaiser a and W. Weise a ' 6 



o 
o 

(N 



oo 



a Physik Department, Technische Universitat Miinchen, D-85747 Garching, Germany 
b ECT*, 1-38050 Villazzano (Trento), Italy 



Abstract 



Qh The role of chiral (pion) dynamics in nuclear matter is reviewed. Contributions to the 



energy per particle from one- and two-pion exchange are calculated systematically, and it 
is demonstrated that already at order in the Fermi momentum, two-pion exchange pro- 
duces realistic nuclear binding together with very reasonable values for the compressibility 
and the asymmetry energy. Further implications of these results are discussed. 



> 

in 
o 

2 1 Introduction 

o t 

\ The present status of the nuclear matter problem is that a quantitatively successful description 
can be achieved, using advanced many-body techniques |], in a non-relativistic framework 
Zj ' when invoking an adjustable three-body force. Alternative relativistic mean field approaches, 
including non-linear terms with adjustable parameters, are also widely used for the calculation 
of nuclear matter properties and finite nuclei ||. At a more basic level, the Dirac-Brueckner 
method || solves a relativistically improved Bethe-Goldstone equation with one-boson exchange 
^ ' AA-interactions. 

In recent years a novel approach to the AA-interaction based in effective field theory (in 
particular, chiral perturbation theory) has emerged j|, [5| . The key element is a power counting 
scheme which separates long- and short- distance dynamics. Methods of effective field theory 
have also been applied to systems of finite density || . 

The purpose of this presentation is to point out the importance of explicit pion dynamics 
in the nuclear many-body problem. While pion exchange processes are well established as 
generators of the long and intermediate range AA-interaction, their role in nuclear matter is 
less evident. The one-pion exchange Hartree term vanishes identically, and the leading Fock 
exchange term is small. Two-pion exchange mechanisms are commonly hidden behind a purely 
phenomenological scalar ("sigma"-) mean field which is fitted to empirical data but has no 
basic justification. This is an unsatisfactory situation which calls for a deeper understanding. 
We report on steps and thoughts in this direction, following ref. [7J. Our approach is closely 
related to the work of Lutz et al. in ref. [HI. 
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Before passing on to our calculation it is useful to draw attention to the following fact. 
A simple but realistic parametrization of the energy per particle, E(kf) = E/A, of isospin 
symmetric nuclear matter is given in powers of the Fermi momentum ^ as 

where the nucleon density is p = 2k^/3n 2 as usual, and M = 0.939 GeV is the free nucleon mass. 
The first term is the kinetic energy of a Fermi gas. Adjusting the (dimensionless) parameters 
a and (3 to the equilibrium density, po = 0.16 fm~ 3 (kfo = 1-33 fm~ l ) and E = E(kf ) = 
— 16MeV, gives a = 5.27 and /3 = 12.22. The compression modulus K = kj (d 2 E(kf)/dkj) kfo 
is then predicted at K = 236 MeV, well in line with empirically deduced values, and the density 
dependence of E(kf) using eq. (1) is remarkably close to the one resulting from the realistic 
many-body calculations of the Urbana group |§. 



2 Chiral in-medium perturbation theory 

The tool to investigate the implications of spontaneous and explicit chiral symmetry breaking 
in QCD is chiral perturbation theory. Observables are calculated within the framework of 
an effective field theory of Goldstone bosons (pions) interacting with the lowest-mass baryons 
(nucleons). The diagrammatic expansion of this low-energy theory in the number of loops has a 
one-to-one correspondence to a systematic expansion of observables in small external momenta 
and the pion (or quark) mass. 

In nuclear matter, the relevant momentum scale is the Fermi momentum kf. At the empirical 
saturation point, &/o — 2m ni so the Fermi momentum and the pion mass are of comparable 
magnitude at the densities of interest. This immediately implies that pions must be included 
as explicit degrees of freedom: their propagation in matter is relevant. Pionic effects cannot be 
accounted for simply by adjusting coefficients of local NN contact interactions. 

Both kf and are small compared to the characteristic chiral scale, 4:irf n ~ 1.2 GeV, 
which involves the pion decay constant /„- = 0.092 GeV. Consequently, the equation of state 
of nuclear matter as given by chiral perturbation theory will be represented as an expansion in 
powers of the Fermi momentum. The expansion coefficients are non-trivial functions of kf/m^, 
the dimensionless ratio of the two relevant scales inherent to the problem. 

The chiral effective Lagrangian generates the basic pion-nucleon coupling terms: the Tomo- 
zawa- Weinberg irnNN contact vertex, (l/4/|)(g^ — q^jl^abc^c, and the pseudovector ttNN 
vertex, [g aI^J -k)^! vH^ai where q a ,b denotes (outgoing) pion four-momenta and qa is the axial 
vector coupling constant (we choose qa = 1-3 so that the Goldberger-Treiman relation = 
gAM/ f n gives the empirical ttN coupling constant, g n N = 13.2). 

The only new ingredient in performing calculations at finite density (as compared to evalua- 
tions of scattering processes in vacuum) is the in-medium nucleon propagator. For a relativistic 
nucleon with four-momentum = (po,p) it reads 

U+M) [ p2 _M2 + l£ - 2 ^P 2 ~ M2 )°(Ponk f - \p\)\ . (2) 

The second term is the medium insertion which accounts for the fact that the ground state of 
the system has changed from an "empty" vacuum to a filled Fermi sea of nucleons. Diagrams 
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can then be organized systematically in the number of medium insertions, and an expansion is 
performed in leading inverse powers of the nucleon mass, consistently with the fcj-expansion. 

Our "inward-bound" strategy [[7] is now as follows. One starts at large distances (small kf) 
and systematically generates the pion-induced correlations between nucleons as they develop 
with decreasing distance (increasing kf). The present calculations are performed to 3-loop 
order (including terms up to order kf) and incorporate one- and two-pion exchange processes. 
The procedure involves one single momentum space cutoff A which encodes dynamics at short 
distances not resolved explicitly in the effective low-energy theory. This cutoff scale A is the 
only free parameter which has to be fine-tuned. (Alternatively, and equivalently, one could 
use dimensional regularization and introduce short-distance physics through adjustable NN 
contact terms). 

We now outline the leading contributions to the energy per particle E{kf). The kinetic 
energy including first order relativistic corrections is 



Ekin{kf) 



10M 



1 - 



hk) ' 
28M 2 



(3) 



Terms of order kf are already negligibly small. At least from this perspective, nuclear matter 
is a non-relativistic system. 

Nuclear chiral dynamics up to three- loop order introduces the diagrams, Fig. 1. They 
include the one-pion exchange (OPE) Fock term, iterated OPE and irreducible two-pion ex- 
change. 





Fig. 1: In-medium chiral perturbation theory: One-pion exchange Fock term (upper left), 
iterated one-pion exchange (upper middle and right) and examples of irreducible two-pion ex- 
change terms. See ref. jffl for details. 

Medium insertions are systematically applied on all nucleon propagators, and the relevant 
loop integrations yield results which can be written in analytic form for all pieces. 
The OPE Fock term becomes 



E lw (k f ) 



9a™1 



(4) 
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where F and G are functions of the dimensionless variable kf/m^. They are given explicitly in 
ref. J7|. All finite parts of iterated OPE and irreducible two-pion exchange are of the generic 
form 



E 2n (k f ) 



mz 



(4*7,0 



k 



k 



g\MH A 1^-1+ m w H 5 [ ^~ 



(5) 



with the functions again given explicitly in ref. |7j . All power divergences specific to cutoff 
regularization are summarized in the expression 



E A (k f ) 



Akj 



■10g A M+(3g A +l){g 2 A -l)A], 



(6) 



(4ttA) 41 

where the attractive and dominant first term in the brackets arises from iterated OPE. Note 
that this term could have been generated, equivalently, by a iViV contact interaction with 
appropriate coupling strength. 



3 Results 

3.1 Nuclear matter equation of state 

A striking feature of the chiral dynamics approach is the simplicity of the saturation mechanism 
for isospin-symmetric nuclear matter. Before turning to the presentation of detailed results, it 
is instructive first to discuss the situation in the exact chiral limit, m n = 0. The basic satu- 
ration mechanism can already be demonstrated by truncating the one- and two-pion exchange 
diagrams at order kj. We can make straightforward contact with the parametrization (1) of the 
energy per particle and identify the coefficients a and (3 of the kj and kj terms, respectively. 
The result for a in the chiral limit is: 

a = ±^(!b!l)*-(!l!l)\ ( 7) 



M \ 4vr J V 4vr 

where we have neglected the small correction proportional to A 2 in eq. (6). The strongly 
attractive leading term in eq. (7) is accompanied by the (weakly repulsive) one-pion exchange 
Fock term. 

The /^-contribution to E(kf) would lead to collapse of the many-body system. The stabi- 
lizing kj-term is controlled by the coefficient (calculated again in the chiral limit) 

p = 4 4 (47r2 + 237 " 24 ln 2) " h = 13 - 55, (8) 

a unique and parameterfree result to this order. Here the two-pion exchange dynamics produces 
repulsion of just the right magnitude to achieve saturation: the result, eq. (8), is within 10% 
of the empirical (3 = 12.2. Adjustment of the short-distance scale A between 0.5 and 0.6 GeV 
easily leads to a stable minimum of E(kf) in the proper range of density and binding energy. 

The full 3-loop chiral dynamics result for E{kf) in symmetric nuclear matter, using m n = 
135 Me V (the neutral pion mass), is shown in Fig. 2 together with a realistic many-body 
calculation. The outcome is remarkable: with one single parameter A = 0.65 GeV fixed to the 
value Eq = — 15.3 MeV at equilibrium, perturbative pion dynamics alone produces an equation 
of state which follows that of much more sophisticated calculations up to about three times the 
density of nuclear matter. The predicted compression modulus is K = 255 MeV, well in line 



with the "empirical" K = (250 ± 25) MeV deduced in refs. |§ |0 . 
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Fig. 2: Energy per particle, E{kf), of symmetric nuclear matter derived from chiral one- 
and two-pion exchange (solid line) j^. The cutoff scale is A = 646 MeV. The dashed line is 
the result of ref. /f|/. 






Fig. 3: Nuclear matter saturation point (E , p ) at finite pion mass = 135 MeV (solid 
line) and in the chiral limit = (dashed line) as function of the cutoff scale A (given 
in MeV). The inserted rectangle corresponds to the empirical saturation point including its 
uncertainties. 
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It is interesting to examine the variation with the pion mass in this context. Fig. 2 displays 
the position of the nuclear matter saturation point (E ,p ), first in the chiral limit and then 
using the physical m n , along lines with varying short- distance scale parameter A. Evidently, 
explicit chiral symmetry breaking by m n is not a qualitatively decisive feature for saturation; 
it influences, however, the quantitative fine-tuning of A. 



3.2 Asymmetry energy 

The specific isospin dependence of two-pion exchange should have its distinct influence on the 
behavior of asymmetric nuclear matter, with increasing excess of neutrons over protons. We 
introduce as usual the asymmetry parameter 5 = (p n — p p )/p = (N — Z)/(N + Z), keeping 
the total density p = p n + p p = 2k^/3n 2 constant. The proton and neutron densities are 
Pp,n — kp n /37r 2 in terms of the corresponding Fermi momenta. Without change of any input, 
we have calculated [0] the asymmetry energy A{kj) defined by 

E as (k p ,k n ) = E{kf) + 5 2 A{k f ) + ... (9) 

The result at nuclear matter density is A = A(kf ) = 33.8 MeV. This is in very good 
agreement with the empirical value A = 33.2 MeV derived from extensive fits to nuclide 
masses [[□] . 

Extrapolations to higher density work roughly up to p ~ 1.5 po- At still higher densities, 
there are indications that non-trivial isospin dependence beyond one- and two-pion exchange 
starts to play a role. A similar statement holds for pure neutron matter which is properly 
unbound, but its predicted equation of state starts to deviate from that of realistic many-body 
calculations at neutron densities larger than 0.2 fm~ 3 . 



3.3 Nuclear mean field from chiral dynamics 

The in-medium three-loop calculation of the energy per particle defines the (momentum de- 
pendent) self-energy of a single nucleon in nuclear matter up to two- loop order. The real part 
of the resulting single particle potential in isospin-symmetric matter at the saturation point, 



for a nucleon with zero momentum, comes out as 12 



U(p = 0, k f0 ) = -53.2 MeV, (10) 

using exactly the same one- and two-pion exchange input that has led to the solid curve in Fig. 
2. The momentum dependence of U(p, kfo) can be rephrased in terms of an average effective 
nucleon mass M* ~ 0.8 M at nuclear matter density, and the imaginary part of the potential 
for a nucleon-hole at the bottom of the Fermi sea is predicted to be about 30 MeV. All these 
numbers are remarkably close to the empirically deduced ones. 



4 Summary and outlook 

Explicit pion dynamics originating from the spontaneously broken chiral symmetry of QCD 
is an important aspect of the nuclear many-body problem. In-medium chiral perturbation 
theory, with one single cutoff scale A ~ 0.65 GeV introduced to regularize the few divergent 
parts associated with two-pion exchange, gives realistic binding and saturation of nuclear matter 
already at three-loop order. At the same time it gives very good values for the compression 
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modulus and the asymmetry energy. These are non-trivial observations, considering that it all 
works with only one adjustable parameter which encodes unresolved short-distance dynamics. 
Of course, questions about systematic convergence of the in-medium chiral loop expansion still 
remain and need to be explored. 

In view of the relevant scales in nuclear matter, the importance of explicit pion degrees of 
freedom does not at all come unexpected. Many of the existing models ignore pions, however. 
They must introduce purely phenomenological scalar fields with non-linear couplings and freely 
adjustable parameters in order to simulate two-pion exchange effects. 

Finally, in order to discuss possible contacts with relativistic nuclear mean field phenomenol- 
ogy, the following working hypothesis suggests itself as a guide for further steps. Assume that 
the nuclear matter ground state represents a "shifted" QCD vacuum characterized by strong 
vector (V) and scalar (S) condensate fields acting on the nucleons, with V ~ —S ~ 0.3 GeV as 
suggested e. g. by in-medium QCD sum rules ]13[. Such a scenario would not produce binding 
all by itself, but establish strong spin-orbit splitting and approximate pseudo-spin symmetry 
14| . Binding and saturation would then result from the pionic (chiral) fluctuations around this 



new vacuum. 



References 

[1] Akmal, A., Pandharipande, V.R. and Ravenhall, D.G., Phys. Rev. C 58 (1998) 1804, and 
refs. therein. 

[2] Ring, P., Prog. Part. Nucl. Phys. 37 (1996) 193, and refs. therein. 

[3] Brockmann, R. and Machleidt, R., Phys. Rev. C 42 (1990) 1965, and refs. therein. 

[4] Kaplan, D.B., Savage, M.J. and Wise, M.B., Nucl. Phys. B 534 (1998) 329; Epelbaum, 

E., Glockle, W. and Meifiner, U.-G., Nucl. Phys. A 671 (2000) 295; Lutz, M., Nucl. Phys. 

A 677 (2000) 241. 

[5] Kaiser, N., Brockmann, R. and Weise, W., Nucl. Phys. A 625 (1997) 758; Kaiser, N., 
Gerstendorfer, S. and Weise, W., Nucl. Phys. A 637 (1998) 395. 

[6] Furnstahl, R.J., Steele, J.V. and Tirfessa, N., Nucl. Phys. A 671 (2000) 396; Furnstahl, 
R.J. and Hammer, H.W., Nucl. Phys. A 678 (2000) 277; Lutz, M., Friman, B. and Appel, 
C, Phys. Lett. B 474 (2000) 7. 

[7] Kaiser, N., Fritsch, S. and Weise, W., Nucl. Phys. A (2001), in print, |nucl-th/0105057 . 



[8] Friedman, B., and Pandharipande, V.R., Nucl. Phys. A 361 (1981) 502. 

[9] Blaizot, J.P, Phys. Reports 64 (1980) 171. 
[10] Vretenar, D., et a!., Nucl. Phys. A 621 (1997) 853; Ring, P., private communication. 
[11] Seeger, PA., and Howard, W.M., Nucl. Phys. A 238 (1975) 491. 

[12] Kaiser, N., Fritsch, S. and Weise, W., preprint (2001), |nucl-th/01080T0| , submitted to Nucl. 
Phys. A. 

[13] Cohen, T.D., Furnstahl, R.J., and Griegel, D.K., Phys. Rev. Lett. 67 (1991), 961; Jin, X., 
et al, Phys. Rev. C 49 (1994) 464. 

[14] Ginocchio, J., these proceedings. 



7 



